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Abstract 

A recently constructed model for low lying excitations in bilayer graphene 
exhibits mid-gap, zero energy modes in its Dirac-like spectrum, when a scalar 
order parameter takes a vortex profile. We show that these modes persist 
when the dynamics is extended by a gauge field interaction, which also renders 
finite the vortex energy. The effect of the gauge field on the zero energy wave 
function is to shift the phase of the (damped) oscillatory component of the 
wave function in the absence of the gauge field. 

The old subject of Dirac zero modes and fractional charge [1] revived recently 
owing to the emergence of graphene as an experimentally realized planar substance 
[2] , whose low-energy excitations can be described by a Dirac equation in two spatial 
dimensions [3] . If the material exhibits various dimerization patterns, the effective 
Dirac fields also interact with a homogenous scalar field (order parameter), and 
this gives rise to a gap in the Dirac spectrum. When the scalar field acquires a 
topologically interesting profile, e.g. a vortex, a zero energy, mid-gap state can occur 
with fractional (fermion) charge [1]. 

An early instance of planar Dirac zero modes was found in [1], but no actual 
experimental setting was given. Today graphene, and graphene-like substances, offer 
the possibility of a physical realization. 
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Monolayer graphene consists of a hexagonal, honeycomb lattice, which may be 
presented as a superposition of two triangular sublattices, A and B. In the tight- 
binding approximation, there are two Dirac points. If a particular dimenization 
— called Kekule distortion — occurs, the effective Dirac Hamiltonian also possesses 
an interaction with a scalar field ip |[5]. 
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The 4x4 Dirac Hamiltonian h^ acts on a 4-spinor 
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where (±) refer to the two Dirac points and {A, B) label the sublattices. The vectorial 
quantities p, a and cr are 2-dimensional. The kinetic term ex ■ p does not mix the 
Dirac points; mixing arises through as a consequence of the Kekule distortion. 
Homogenous (p = m produces a mass gap, while an n-vortex profile for (fi{r) produces 
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zero modes. The Hamiltonian ([T]) anti-commutes with 



a^hia^ = ~hi (3) 

Therefore maps positive energy solutions onto negative energy solutions and zero 
modes can be chosen to be eigenmodes of a^. This "energy reflection symmetry" is a 
manifestation of the sublattice symmetry found in the honeycomb graphene lattice. 

The above model was extended by including an interaction with a gauge field. A, 
whose purpose is to unpin the vortices; in the absence of the gauge field the scalar 
vortex carries infinite energy (per unit legnth) [6]. 
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The gauged model, possesses a local chiral gauge symmetry 

A ^ A + Vw, 
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and one readily verifies the identity 



ht = exp (a^ bj hi exp (a^^^b], 
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Thus the extended model still retains the energy reflection symmetry, and possess 
zero energy eigenmodes, whose wave functions differ from those with just a scalar 
vortex by the factor e~" v^" ^ , 

However, it may be difficult to achieve experimentally the Kekule distortion. 
Recently a model that is physically different but mathematically similar to ([T]) has 
been put forward, with the suggestion that the excitation condensate, needed for 
topological effects, fractional charge etc., can "be produced in the laboratory in the 
near future" |7j. The physical system consist of a graphene bilayer, separated by 
a dielectic barrier, and biased by an external, constant voltage V . In a mean-fleld 
approximation, the Hamiltonian for the above bilayer system is given by 
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which acts on the 4-spinor \E'. 
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As before (A, B) refer to the sublattices, but (P, [2]) label the two layers, which are 
nested, one directly above the other. There are no Dirac point labels, because the 
above description refers to a single Dirac point in each lattice of the two stacked 
lattices. Here ip describes the condensate arising from states bound by interlayer 
Coulomb forces between particles in one layer and holes in the other. This dynamics 
is modeled by a 4-Fermi interaction of strenght \J . A gap equation is solved in the 
Hartee-Fock approximation, leading to 
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Here t is the hopping amplitude between sites on each of the two monolayers; there 
is no interlayer hopping within this model's approximations. Eq. f llUI) holds in the 
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limit A ^ V ^ \(p\, where A is an energy cut off It is striking that this order 
parameter enters the bilayer theory in a way identical to the Kekule distortion of the 
monolayer model, /ii in ([T]). 

The presence of 'j^V in h2, which has no analog in hi, spoils the energy reflection 
symmetry ([3]). But another property of /12 ensures similar behavior. One verifies 
that /12 satisfies 

pa^hlf]a^ = h2. (11) 

Thus energy reflection works as 

^.E = f3a^^*E, (12) 
and /12 possesses zero-energy eigenstates, satisfying [7j 

«2/?^o = ^*. (13) 

In view of our earlier work on gauging the monolayer graphene model [B], we are 
led to study the gauged version of /i2. 

/i^ = a • (p - 75 A) + /5 1(^1 e-'^' ^ + 75 1^ (14) 

Gauge transformations follow ([5]) and V is gauge invariant. The new energy reflection 
property, ffTTl) . f|T2l) . is maintained. Consequently we expect to find zero modes, which 
we now exhibit . 

The 4-spinor is presented in terms of 2-spinors. 




With our Dirac matrices, the zero energy spinors satisfy according to (IT5]1 . \l/2 = 
cr^ and the eigenvalue equation reads 

(o--(p-A) + l^) + = (16a) 

^* ^1 - (o- ■ (p - A) + V) = (16b) 

In fact, the second equation is a consequence of the first, and needs not be con- 
sidered separately. Continuing with the matrix reduction, we set ^I/i(r) = 
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and eq. (]16ap now reads 
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where we have taken = —ne^^ ^ A (r), with ^4(0) = 0,A(oo) = | and ip = 
m (r) e*"^, with m (0) = 0, m(oo) = m. 

To separate the angular dependence, to make the equations real and to simphfy 
them, we posit the Ansatz 
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/i = I2 = M'(r) = m(r) and /, g are real. Single valuedness requires that n 
be an odd integer. The final equations read 
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When A remains unspecified (apart from its asymplotes) eqs. (fT9!l do not appear 
explicitly integrable. Nor can A be removed, as in the monolayer case (Q, ([7]). 
However, one can show that a normalizable solution exists. 

For r — i> 00, A — i> I and (HM reduce to 



f -Vg = 0, 

g' + vf = 0, 

with solution that involves two constants, {c,d). 

f{r) = c cos Vr + d sin Vr 
g{r) = —c sin Vr + d cos Vr 
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Evidently owing to the r~5 e"^^*^'') factor both F and G are always damped at large 
r. Thus the wave function will be acceptable and normalizable if a solution that is 
regular at the origin can be constructed. 

At the origin A vanishes and the equations (fl9|) reduce to 

d n 
dr 2r 



f-Vg = (22a) 



| + ^)s + V-/ = (22b) 

Of course these are the same equations, which hold for all r in the absence of A, as 
with the Hamiltonian /12 in ([8]). Their solution is given in terms of Bessel functions 

f = A Jn_i_ iVr), g = -A J„^i {Vr) (23) 

Note that the large r asymptote of (123!) is of the form (I2T!) with specific values for 
c = \ - cos = \ - sin Thus the effect of the gauge field is to move c and 

yvr 4'y7r4 oo 

d from the above values; i.e. A causes a phase shift in the profiles without gauge 
field. 
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